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A crucial step in revealing the nature ol unconventional superconductivity is to investigate the 
symmetry of the superconducting order parameter. Scanning tunneling spectroscopy has proven a 
powerful technique to probe this symmetry by measuring the quasiparticle interference (QPI) which 
sensitively depends on the superconducting pairing mechanism. A particularly well suited material 
to apply this technique is the stoichiometric superconductor LiFeAs as it features clean, charge 
neutral cleaved surfaces without surface states and a relatively high Tc ~ 18 K. Our data reveal 
that in LiFeAs the quasiparticle scattering is governed by a van-Hove singularity at the center of 
the Brillouin zone which is in stark contrast with other pnictide superconductors where nesting is 
crucial for both scattering and s±-superconductivity. Indeed, within a minimal model and using the 
most elementary order parameters, calculations of the QPI suggest a dominating role of the hole- 
like bands for the quasiparticle scattering. Our theoretical findings do not support the elementary 
singlet pairing symmetries s++, s±, and d-wave. This brings to mind that the superconducting 
pairing mechanism in LiFeAs is based on an unusual pairing symmetry such as an elementary p- 
wave (which provides optimal agreement between the experimental data and QPI simulations) or a 
more complex order parameter (e.g. s -\- id- wave symmetry). 



The discovery of iron-based superconductors [T] has 
generated enormous research activities to reveal the na- 
ture of superconductivity in these materials. In particu- 
lar, s±-pairing, i.e. an s-wave order parameter with al- 
ternating sign between almost perfectly nested hole and 
electron pockets has been suggested to be prevailing for 
the entire class of iron-based superconductors fS'-^. In 
this regard, the material LiFeAs is of particular inter- 
est since experiments have proven an absence of nesting 
[S^ and theoretical works yield contradictory results, i.e., 
both s-t-wave singlet as well as p-wave triplet pairing has 
been suggested O [7] . 

A powerful method to probe the symmetry of the su- 
perconducting order parameter is to map out the spatial 
dependence of the local density of states (DOS) by scan- 
ning tunneling spectroscopy (STS). In such experiments, 
the relation D0S(i5) cx d//(iV(Vbias) with tunneling volt- 
age V and current / at energy E and bias voltage Vbias 
{E = eVbias) is exploited. The spatial dependence of the 
DOS often arises from an impurity scattering of the con- 
duction electrons. In this case the incident and scattered 
quasiparticle waves interfere and give rise to Friedel-like 
oscillations in the local density of states (LDOS). Such 
quasiparticle interference (QPI) effects have first been 
observed by STS experiments on normal state metal sur- 
faces [8l-fl0] . A convenient way to extract the dominating 
scattering vectors q from a spatially resolved image of 
the QPI pattern is by analysis of its Fourier transformed 
image [TU] . This so-called spectroscopic-imaging scan- 
ning tunneling microscopy (Sl-STM) has proven to be an 
extremely powerful technique to investigate unconven- 
tional superconductors because the scattering rate be- 
tween quasiparticle states with momenta k and k' is pro- 



portional to coherence factors (wkw^, =F VkV^, ) , where the 
=p sign is determined by the nature of the underlying scat- 
tering mechanism. The coherence factors are sensitive to 
the phase of the superconducting order parameter via 
the Bogoliubov coefficients Wk and Wk which fulfill the 
relation fk/wk = {E]^ — ^k)/A^ with the quasiparticle 
energy Ek = ±(^^ + \Ak\Y^^ jUllIl]. Thus, through 
the coherence factors the QPI pattern is decisively influ- 
enced by the nature of superconductivity, in particular 
the symmetry of the superconducting gap. Prominent 
examples are recent studies on the QPI of cuprate high- 
temperature superconductors |13H15j and, more recently, 
also of iron-based superconductors [IS] . 

LiFeAs is a stoichiometric superconductor which ex- 
hibits clean, charge neutral cleaved surfaces with a bulk- 
like electronic structure at the surface [17]. Its critical 
temperature is Tc ~ 18 K [18] and currently available 
single crystals are of high quality [TO] . With these prop- 
erties LiFeAs emerges as a prime material to probe the 
superconducting gap characteristics using SI-STM. Our 
studies reveal that the QPI is governed by a van-Hove 
singularity at the center of the Brillouin zone where the 
scattering between hole-like bands dominates. Moreover, 
the comparison of the experimental data with QPI calcu- 
lations suggests incompatibility with elementary singlet 
pairing symmetries (s±, s++ or d-wave) but yields opti- 
mal agreement if an elementary p-wave order parameter 
is assumed. 

We have grown high quality single crystals of LiFeAs 
by a self-flux technique [19j. The SI-STM measurements 
have been performed in a home-built apparatus at T w 
5.8 K in cryogenic vacuum [20] . 

Figure [T^ shows a representative surface topography 
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FIG. 1: (a) Surface topography of LiFeAs measured in con- 
stant current mode (7 — 600 pA, Vbias ~ —50 mV) taken at 
T ~ 5.8 K. Black arrows indicate the direction of the lattice 
constants [T5] with a — 3.7914 A. (b) Spatially averaged tun- 
neling spectrum taken in the square area (dashed lines) in 
(a). The spectrum exhibits a gap with 2 A ~ 10 mV, which 
(taking thermal broadening into account) is consistent with 
low-temperature (~ 500 mK) tunneling data [21] of the su- 
perconducting gap of LiFeAs. Representative energy values 
for further QPI analysis in Fig. [2] are marked by red circles. 

of LiFeAs after cleaving in cryogenic vacuum at 5.8 K, 
registered in constant-current mode at the bias voltage 
Vbias = —50 mV, evidencing the pristine nature of the Li- 
terminated surface. The field of view exhibits a highly 
periodic atomic corrugation vifhich corresponds to about 
2000 atoms in the topmost layer. 14 defects are present 
in the field of view and appear as valleys superimposed on 
the atomic corrugation. This yields a defect concentra- 
tion < 1% and highlights the high purity of our samples. 

In order to resolve the QPI we have measured 
the spatial- and energy resolved differential tunnel- 
ing conductance, dI/dV{r, Vbias) in the range E — 
—50 meV ... 50 meV in full spectroscopic mode. A spa- 
tially averaged dl/dV spectrum, taken on a defect-free 
area clearly reveals a superconducting gap (Fig. [T]d) . SI- 
STM maps at representative energies (Fig. [2^-h |^) re- 
veal very clear QPI patterns which are most pronounced 
in the vicinity of the coherence peaks at negative en- 
ergy (Vbias ?J —20 mV). In this energy range, the QPI is 
clearly not only visible as relatively strong modulations 
close to the defects but also appears as clear wave-like 
modulations (with a wavelength of a few lattice spacings) 
in the relatively large defect-free area in the center of the 
field of view. QPI patterns are also discernible at positive 
energy, but compared to the pronounced modulations at 
negative values, the amplitude of the modulations decay 
more rapidly when moving away from a defect. 

Figure [2j-p displays the Fourier transformed images 
of the shown conductance maps. The most salient fea- 
ture is a bright structure distributed around q = (0,0). 
In similarity to the observed real-space modulations this 
feature is particularly pronounced at energies Vbias ~ 
[—20 mV, 0] where it attains a squarish shape with the 
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FIG. 2; (a-h) SI-STM maps of the region shown in Fig. [T] 
at selected representative bias voltages, (i-p) Fourier trans- 
formed images of the maps shown in a-h. Bright spots at 
(±7r, ±7r) and at higher q result from the atomic corrugation 
in the real space images. 



corners pointing along the {qx,0) and {0,qy) directions. 
Upon increasing Vbias to positive values, the intensity 
at the square corners increasingly fades and for Vbias > 
10 mV the squarish shape changes to an almost round 
structure which remains in that shape up to 50 mV. The 
Fourier transformed images also reveal further well re- 
solved structures with significantly lower intensity cen- 
tered around (7r/2, 7r/2), (tt, 0) and (vr, n) which again are 
most pronounced between -20 mV and the Fermi level. 
At more negative bias, these finer structures fade, while 
at positive bias voltage they develop into a rather fea- 
tureless diffuse background. 

Already without considering a full calculation of the 
QPI in LiFeAs we can deduce important qualitative in- 
formation from the data. In fact, the observed Fourier 
transformed images possess an eye-catching similarity to 
the Fermi surface of LiFeAs [51 [7], which immediately 
implies that the QPI is dominated by a set of scatter- 
ing vectors which emerge from one relatively small, high 
DOS region in the Brillouin zone. This is illustrated in 
more detail in Figure[3^ and b where we compare the con- 
stant energy contour (CEC) of LiFeAs at E = —11.7 mV 
[5] with the observed QPI intensities in the Fourier trans- 
formed image. Most prominent is that the observed cen- 
tral squarish structure in Figure appears like a some- 
what enlarged smeared replica of the large hole-like CEC 
around (0,0) [5S]. This observation can directly be un- 
derstood as stemming from scattering processes (qi) con- 
necting the very small and the larger squarish-shaped 
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FIG. 3: (a) Simplified GEO [5] at _E = —11.7 meV in the periodic-zone scheme of the Brillouin Zone (BZ), where the first BZ 
(referring to the unit cell with two Fe atoms) is indicated by the dashed lines. However, the used coordinates in reciprocal space 
refer to the unit cell with one Fe atom, in order stay consistent with the theoretical work in Ref. [T] This choice of reciprocal 
coordinates leads to Bragg-intensity at (±7r, zbvr) instead of (±27r, 0) (and (0, ±27r)) as one would expect for a two-Fe unit cell. 
The two pockets around (0,0) represent hole-like CEC while the pockets at the zone boundary are electron-like. qi,2 represent 
scattering processes which connect states on the small hole-like GEG and on other CEG, qs and q4 represent scattering between 
the electron-like and within the large hole-like GEG, respectively, qs^e represent umklapp processes. Note that each scattering 
process qi,...,6 is described by a set of scattering vectors as is illustrated for qi (dashed and solid arrows), (b) Measured Fourier 
transformed image at the same energy (the same as Figure 21) with qi,...,6 superimposed. The most salient QPI features around 
(0,0) and (tt, 0) match well with qi and q2 (see text). The further observed but less prominent QPI intensities around (tt, tt) 
and at (7r/2,7r/2) are well described by qs and q4, respectively. The umklapp scattering vectors qs and qe might also be of 
relevance here, (c-g) Galculated QPI in q space assuming the normal state and a superconducting order parameter with s±-, 
d-, p-, and s++ symmetry. 



hole- like CEC around (0,0). Furthermore, the much 
weaker structure at q = (I'jO) in Figure 3b, which re- 
sembles a smeared version of the electron-like CEC near 
k = (tt, 0), apparently results from scattering processes 
(q2) connecting this CEC with again the small hole-like 
CEC. Note that due to the presence of a van-Hove singu- 
larity [5] at k = (0, 0) the DOS for momentum vectors on 
the small hole-like CEC is enhanced significantly. Thus, 
scattering processes combining those momenta with k- 
points close to the other CEC are dominant and give rise 
to a Fourier transformed QPI image which essentially 
represents a 'map' of all CEC at a given energy. 

Having assigned the main features of the measured QPI 
patterns to the normal state electronic structure, the next 
step is to establish how in the superconducting state the 
DOS is redistributed by the opening of the superconduct- 
ing gap. One expects that in the superconducting state 
the DOS at energy values close to the gap value is fur- 
ther boosted in comparison to the normal state since the 
quasiparticle dispersion £^k = ±(?k + l^kP)^^^ is rather 
flat in the vicinity of the Fermi surface. Furthermore, 
depending on the gap function Ak particular scattering 
channels are suppressed while others are enhanced ac- 
cording to the coherence factors. Consequently, the QPI 
measured at energies \E\ close to the averaged gap value 
will be enhanced significantly as compared to the nor- 
mal state. This intensity redistribution contains detailed 
information about the structure of the superconducting 
order parameter. 

For a quantitative analysis of the experimental data, 
we calculated the QPI in the superconducting state using 
an appropriate BCS model for LiFeAs which can describe 
three cases of elementary singlet pairing (s++, s±, and 
d-wave) as well as a p-wave triplet pairing scenario. To 
ensure that our calculations are based on unbiased exper- 



imental findings, the electronic band structure of LiFeAs 
was modeled using tight-binding fits derived from recent 
ARPES experiments [5]. Wc note that our band struc- 
ture model obtained in this way is, at least for the one- 
particle states to be considered here, consistent with an 
alternative three band orbital model [7] for LiFeAs. Fol- 
lowing Refs. and [23J we numerically evaluated the 
LDOS in momentum space. The scattering potential 
of the impurities was modeled by a conventional non- 
magnetic local Coulomb potential PO] . 

In Fig. [3j:-g we compare the results of our calcula- 
tions for the energy E = —11.7 niV. The calculated QPI 
pattern of the normal state shown in Figure |3j: exhibits 
an extended low-intensity region around (0,0), in stark 
contrast to the experimental data (Fig. ^jp). The situ- 
ation drastically changes if the superconducting state is 
taken into account in the QPI calculations. Figure [3ji- 
f shows the calculated QPI for four different types of 
the superconducting order parameter. The main feature, 
which now can be observed for three of the four supercon- 
ducting cases is a wide high-intensity region centered at 
(0,0). This qualitative change clearly corroborates that 
the intense central structure observed in the experiment 
is characteristic of the superconducting state and a com- 
parison of Figure [3^ with the experimental result clearly 
allows us to exclude a clean d-wave pairing in LiFeAs. 
Furthermore, from the location of the high intensity re- 
gion we conclude that in the superconducting state the 
dominating scattering processes are those which connect 
states on the hole-like bands. The prevalence of these 
scattering processes is naturally explained by the pres- 
ence of the van Hove singularity in the center of the Bril- 
louin zone and is one of our major findings. It implies 
that nesting-enhanced scattering, which is predominant 
and crucial for s± -superconductivity in other iron pnic- 
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tides [im [ini [22] , is unimportant in LiFeAs. 

A closer examination of the individual QPI patterns 
brings to light a very strong dependency of the pattern 
on the assumed pairing symmetry. This allows a direct 
comparison with the experimental result. Interestingly, 
we find a striking agreement between the experimentally 
observed QPI and the calculated image of Fig. |3f, which 
is based on a triplet paired superconducting state. In 
this case we considered a gap function which has a chiral 
p-wave symmetry to be consistent with the experimental 
observation of lacking gap nodes 0. The agreement is 
obviously much less pronounced for s±-, S-|_+ and d-wave 
singlet pairing. The main bright feature around (0, 0) 
is absent for the d-wave or rather circular for the s_(_+- 
wave cases where we have included higher harmonics in 
order to account for the experimentally observed varia- 
tion of the gap size along the Fermi surface. One could 
argue, that the intensity distribution around (0, 0) found 
for the assumed s-t -pairing somewhat resembles the ex- 
perimental ndings. However, the well resolved structures 
which are described by the scattering vectors q2 and qs^e 
which are marked by green and yellow lines in Fig. |4f^a) 
are consistent with the experimental result only for the 
p-wave case. Thus, among the three considered order pa- 
rameters the best agreement between experimental data 
and calculated QPI images is apparently obtained for a 
p-wave order parameter. 

To show that this surprising similarity holds also at 
other energies, we present in Fig. |4] a comparison be- 
tween experimental results and calculated p-wave based 
QPI for another four energy values. We find an excellent 
agreement for the energy evolution of the three most in- 
tensive QPI structures (framed in Fig. |4]ja)) which are 
related to the scattering vectors qi, q2 and qs^g (com- 
pare Fig. [Sjja)). A closer inspection of the boundary of 
the central structure which is shown in Fig. |4](b) reveals 
a clear dispersive behavior of the high-intensity region 
within some finite energy range. As it is presented in 
Fig. gc) this energy dependence is confirmed very well 
within our model calculations. 

Our main observation that the quasiparticle scattering 
is dominated by the van-Hove singularity at the center 
of the Brillouin zone fits in very well in recent exper- 
imental and theoretical findings. Angle resolved photo 
emission spectroscopy (ARPES) data provide clear-cut 
evidence for the absence of strong nesting in LiFeAs [5] . 
From these experiments an s± pairing scenario, which is 
believed to be nesting- driven by antiferromagnetic fluctu- 
ations at q = (tt, 0) 0, is not expected. In fact, a recent 
theoretical analysis based on a realistic band structure 
model with poor nesting yields a dominating role of the 
small hole-like Fermi surface for the scattering in LiFeAs, 
giving rise to ferromagnetic fluctuations which drive an 
instability towards triplet superconductivity . It seems 
rather encouraging that our QPI calculations which are 
based on a quite simple model and very elementary as- 
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FIG. 4: (a) Comparison of the measured QPI pattern and 
its calculated counterpart based on a chiral p-wave supercon- 
ducting order parameter vor various energies. The clearly ob- 
servable dispersion unambiguously demonstrates that the ob- 
served structures result from QPI and not from non-dispersing 
modifications of the DOS around a defect, b) Representative 
experimental scans along q = {h, h) with Gaussian fits which 
yield the edge position of the measured central structure along 
this direction, c) Size of the measured and calculated central 
structure along q = {h, h). 



sumptions (potential scattering, i.e., no specific presump- 
tions for the scattering centers as well as elementary or- 
der parameters) are consistent with this finding. It is, 
however, thinkable that a more complex order parame- 
ter (such as s + id symmetry) and/or a specific character 
of the impurities would yield a similar convincing agree- 
ment. Thus, further experimental and theoretical studies 
are indispensable for clarifying the nature of supercon- 
ductivity in LiFeAs. 
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Experimental details 



Model of the Quasiparticle Interference 



High quality single crystals of LiFeAs with a critical 
temperature Tc ~ 18 K have been grown by a self-flux 
technique |19| . The chemical and physical properties of 
the crystals were studied by means of energy-dispersive 
X-ray (EDX) analysis, inductively coupled plasma mass 
spectroscopy (ICP MS), single crystal diffractometry, nu- 
clear quadrupole resonance (NQR), magnetization, spe- 
cific heat, electrical transport and angle resolved photo 
emission spectroscopy (ARPES) [51 [TO], 

The crystals were inserted into a home-built low- 
temperature scanning tunneling microscope and cooled 
down to T « 5.8 K where the samples were cleaved 
in cryogenic vacuum. Subsequently scanning tunneling 
microscopy (STM) and -spectroscopy (STS) measure- 
ments were performed. After cleavage the sample sur- 
face exhibited extended areas with atomic corrugation. 
A rather high density of further objects, presumably left- 
over atoms from the cleaving process, were present in the 
first topographic measurement after cleavage. A part of 
these objects were highly mobile and were 'pushed' out 
of the investigated area by extensive scanning. This re- 
sulted in the topographic data shown in Fig. 1, which 
displays an atomically flat surface containing the immo- 
bile impurities which give rise to the quasiparticle inter- 
ference. 

The SI-STM maps shown in Fig. 2 have been acquired 
in full spectroscopic mode, i.e., by measuring STS at each 
of the 256 x 256 pixel by stabilizing the tip with feedback 
loop on at a setpoint of Vbias = — 50 mV and / = 600 pA 
followed by subsequently ramping Vbias to -1-50 mV with 
the feedback loop switched off. During ramping the volt- 
age we recorded /(14ias) and dl/dV (Vbias) where the lat- 
ter was recorded using a lock-in amplifier with a modu- 
lation voltage Knod = 1-2 mV (RMS) and a modulation 
frequency /mod = 3.333 kHz. 

The Fourier transforms in Fig. 2i-p were calculated 
from the SI-STM maps in Fig. 2 a-h after filtering the 
data using a Hamming window in order to suppress edge 
effects. Since the resulted Fourier transformed images 
were slightly distorted due to piezo relaxation effects, the 
images were first linearly transformed in order to restore 
the tetragonal lattice. Afterwards each Fourier trans- 
formed image was symmetrized with respect to the high 
symmetry in-plane directions. 



Following the ideas of Refs. [231 and [23| we developed 
a theoretical model to calculate the tunneling conduc- 
tance which is measured by an STM experiment. Within 
a simplified model of the tunneling process the differen- 
tial tunneling conductance dI/dV{r, V) at a location r 
and Voltage V is proportional to the single-particle den- 
sity of states at energy w. The spatial dependence of the 
tunneling-matrix element, which includes contributions 
of the sample wave function around the tip, has been 
neglected within the calculations. To study the quasi- 
particle interference we calculated the fluctuations of the 
density of states (LDOS) arising from an electron scat- 
tering off impurities which are distributed in the sample. 

The fluctuations of the LDOS are determined by the 
Green's function of a superconductor in the presence of 
impurities. Assuming the electron field inside a super- 
conductor can be described by a four-component vector 
which is usually written in Nambu notation as 

vl't(r,r)= (V'|(r,T),V|(r,T),V'i(r,T),-V'T(r,r)) , (1) 

where r denotes the position vector and r the imaginary 
time, the matrix Green's function is defined as the time- 
ordered average, 

G„^(r', r; r) = -(r,vE'„(r', r)*^ (r, 0)). (2) 

As it is usual for an electronic system in the presence of 
impurities we calculate the Green's function [2] using the 
t-matrix method. For the dominant effect the scatter- 
ing potential is taken into account in Born approxima- 
tion. Within this approximation the Fourier transformed 
Green's function can be written as 

G(k,k',w) = Go(k,w) + Goik,u;)U{k,k',uj)Goik',uj), 

(3) 

where U{k,k') is the scattering potential of a sin- 
gle impurity and Go(k, w) is the bare Green's func- 
tion of the BCS superconductor without impurities. 
The LDOS is determined by the analytic continuation 
G(r', r; za;„) G(r', r; z) of the Matsubara Green's func- 
tion G(r', r; zw„) = G(r', r; r)e*(2n+i)^Tr^^^ 



p{r,uj) = -ImTr^-4r-^ 

TT Z 



G(r, r; w — iS) 



(4) 



where 1 is the 4x4 unit matrix and is the isospin 
matrix. 



T3 



1 Q 

-1 



(5) 



where an underscore denotes a two dimensional matrix. 
Taking the Eqs. [3] and [4] we evaluated the LDOS using 
a BCS model for a triplet superconductor which is based 
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on a realistic three band model for LiFeAs taken from 
Ref. [21 Note that the band structure of the used model 
is consistent with recent ARPES experiments. For the 
superconducting order parameter we assume a chiral p- 
wave symmetry with the following d-vector, 

d(k) = A[sin(fc2;) ± i sin(fcj^)]ez, (6) 

where A is chosen such that the corresponding gap size 
in the calculated conductance dI/dV{uj) coincides with 
the experimental one shown in Fig. |4j Note that within 
our simplified model the conductance is assumed to be 
proportional to the local density of states which is given 



byEq.|4} 

By assuming that the scattering potential of the 
impurities is screened at length scales comparable to 
the lattice spacing we model the impurity scattering 
by a local Coulomb potential which can be described 
by a matrix of the form U = Uf^. Note that for a 
local potential the momentum dependence of [/ can be 
omitted. Furthermore, for simplicity, we assumed that 
the scattering is equal for all electron bands. 



